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We apply pseudo-spectral methods to construct global solutions of functional renormalisation 
group equations in field space to high accuracy. For this, we introduce a basis to resolve both finite 
as well as asymptotic regions of effective potentials. Our approach is benchmarked using the critical 
behaviour of the scalar 0(1) model, providing results for the global fixed point potential as well as 
leading critical exponents and their respective global eigenfunctions. We provide new results for (1) 
multi-critical 0(1) models in fractional dimensions, (2) the three-dimensional Gross-Neveu model 
at both small and large N, and (3) the scalar-tensor model, also in three dimensions. 


I. INTRODUCTION 

Many interesting physical phenomena are charac¬ 
terised by strong coupling. Among them, there are very 
fundamental problems as for example confinement in 
QCD or potentially the quantisation of gravity. Con¬ 
ventional methods from quantum field theory, such as 
perturbation theory, are not applicable in such cases. 
In recent years, an increasingly successful method has 
been the functional renormalisation group (FRG). It is 
based on the Wilsonian idea of integrating out modes 
momentum-shell by momentum-shell. In the last two 
decades, the FRG, based on the formulation by Wet- 
terich [1], was successfully applied to a range of topics, 
including scalar field theories [2-11], fermionic systems 
[12-18], critical phenomena [19-25], gauge theories [26- 
33] and quantum gravity [34-48]. From a technical per¬ 
spective, the generic outcome of applying the FRG to a 
given model is a coupled system of non-linear (integro- 
)differential equations of complex structure. This is due 
to the full field- and momentum-dependent propagator 
entering the Wetterich equation. There exist only a few 
cases where the full equations can be solved analytically. 
In all other cases one has to consider the system within 
some truncation, retaining only a manageable number of 
operators. Even then, the equations are rarely analyti¬ 
cally solvable, e.g. in a large-or mean-field approxima¬ 
tion. However, if one seeks a solution without such ap¬ 
proximations, numerical methods appear indispensable. 
Various methods are used depending on the special struc¬ 
ture of the system of equations. All these methods aim at 
a numerically accurate solution that fulfils the equations 
to high precision. 

In this work, we advocate the use of pseudo-spectral 
methods, especially referring to Ghebyshev polynomials 
as basis functions as a versatile tool for FRG equations. 
The fall-off of the corresponding expansion coefficients 
provides a measure for the accuracy of the approximate 
solution. Pseudo-spectral methods are a well-suited, fast 
means to treat a wide range of different problems: ordi¬ 
nary and partial differential equations as well as eigen¬ 
value problems [49] . These are all problems relevant both 
in the FRG formalism and physics in general. 

Here, we focus on calculating global solutions of phys¬ 
ical systems, especially referring to ordinary differential 


equations (ODEs). In particular, we cover the basics of 
the pseudo-spectral concepts, and ultimately apply them 
to several interesting problems arising from FRG applica¬ 
tions. Some of them, as the Ising model in three dimen¬ 
sions, have already been the focus of many detailed stud¬ 
ies in the past, making them a perfect testing ground for 
the methods presented here. It is convenient to start with 
these models to discuss some mathematical and techni¬ 
cal details of the equations and their implementation. In 
this way, one can also easily compare to known results. 
In order to emphasise that the range of applications of 
the methods presented here is very large, we present a 
variety of applications and extract several new results. 
This includes multi-critical phenomena in non-integer di¬ 
mensions, and the three-dimensional Gross-Neveu model. 
For the scalar-tensor model proposed recently in [50] , we 
gain new insights which cannot be obtained within local 
expansions. Let us also point out that the methods pre¬ 
sented here are heavily used in other contexts [49, 51], 
as e.g. finding solutions to Einsteins equation [52, 53]. 
First applications to FRG problems have been given in 
[54-56]. Our approach is particularly suited for global 
aspects and also resolves asymptotic behaviour in a con¬ 
trolled way. Recently, it has been used to globally resolve 
the supersymmetric analogue of the Wilson-Fisher fixed 
point [57]. 

This work is organised as follows: in section H, we 
present the basic ideas of pseudo-spectral methods, fo¬ 
cusing on the properties essential for the subsequent dis¬ 
cussion. Afterwards, section HI sheds light on the specific 
application of the methods to a given problem, while sec¬ 
tion IV is a short overview on the FRG. We then turn 
our attention to the 0(1) model in section V, first study¬ 
ing the Wilson-Fisher fixed point in three dimensions. 
Furthermore, we extend known results on multi-critical 
fixed points in non-integer dimensions. Gonsequently, 
section VI treats the Gross-Neveu model, first in the large 
flavor number limit, then considering finite flavor num¬ 
bers. Finally, section VH discusses a scalar-tensor model 
which couples a scalar field non-minimally to gravity. 

The numerical results were obtained with code written 
in G-I-+, including the libraries BOOST [58] for handling 
arbitrary precision and Eigen [59] for dealing with linear 
algebra. To analyse and present the data, Mathematica 
10 [60] was employed. 
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II. PSEUDO-SPECTRAL METHODS 

Pseudo-spectral methods aim to represent a function 
via a suitable expansion. Suitable means that the ex¬ 
pansion shall be accurate and easily treatable. In the 
present context, all derivatives needed should be easily 
computable to high accuracy and they and the function 
itself should be evaluable at arbitrary points. Natural 
candidates for such an expansion are the classical orthog¬ 
onal polynomials, which have some convenient properties 
regarding convergence, evaluation and taking derivatives. 
In this paper, we will focus on the Chebyshev polynomi¬ 
als of the first kind, which are defined by 

Tn{cos{x)) = cos{nx), n e No , (I) 

and their cousins, the rational Chebyshev polynomials 

[61], 

Rn{x)=T^(^^y ( 2 ) 

Here, L > 0 is an arbitrary parameter, encoding the 
precise compactification in x. The reason for this choice 
is that they have superior convergence properties as com¬ 
pared to Legendre polynomials or Chebyshev polynomi¬ 
als of the second kind. Both Hermite and Laguerre poly¬ 
nomials are ill-suited for our problems for the following 
reason: they are defined on unbounded intervals, and in¬ 
creasing the interpolation order changes the asymptotic 
behaviour. For the problems usually encountered, the 
asymptotic behaviour is fixed and thus the convergence 
properties of Hermite and Laguerre polynomials is diffi¬ 
cult to control. One possibility is to use the Hermite or 
Laguerre functions, which decay exponentially, but again, 
the convergence properties are often worse compared to 
those of rational Chebyshev polynomials. 

In the following, we will collect some important proper¬ 
ties of Chebyshev polynomials. Similar relations hold for 
their rational counterparts. We will only state results, for 
deeper information and proofs, consider e.g. [49]. First, 
a fast way of evaluating an expansion in Chebyshev poly¬ 
nomials, 

N 

f{x) = '^a^T,{x), (3) 

1=0 

at an arbitrary point x is given by the Clenshaw algo¬ 
rithm, which is based on the 3-term recurrence relation 

Tn+i{x) = 2xTn{x) - Tn-l(x) . (4) 

With this, the function f{x) is readily evaluated via 
the recursive algorithm 

i>N+2 = bjst+i = 0 , 

bi — CLi 2xbi-\-i bi^2 , 

(5) 


Second, the derivative of f(x) can again be expanded 
in a sum of Chebyshev polynomials of degree N —1. The 
expansion coefficients, call them a], are given recursively 
by 

fljv-i = 2Na]s[ , 
ajv-2 = 2(A — l)aAr_i , 

= 2(i -f l)aj+i -b a-+2 • (6) 

Both the Clenshaw and the derivative algorithm are 
high-performance algorithms and numerically stable due 
to their recursive nature. 

The third, but most important property of the Cheby¬ 
shev polynomials is the exceptional convergence of the 
expansion coefficients, which is based to their relation to 
Fourier series as indicated by (I). To be more precise, let 
us first define the algebraic index of convergence as the 
largest number k for which 

lim |a„|n^ < oo . (7) 

n—>-oo 

So-called exponential convergence is achieved if the co¬ 
efficients Qn decrease faster than I/n^ for any k. For 
example, the Chebyshev expansion of a Lipschitz con¬ 
tinuous function is always converging exponentially. One 
can further differentiate exponential convergence into su¬ 
pergeometric, geometric and subgeometric convergence, 
but these details shall not matter here. On a log-log plot, 
algebraic convergence manifests itself in a straight line, 
whereas exponential convergence is indicated by a bend¬ 
ing downwards. Note that the definition of convergence is 
an asymptotic one, and might set in only when including 
a large number of coefficients. 

By Darboux’s principle, convergence properties of a se¬ 
ries is related to the singularity structure of the function 
to be interpolated. This includes poles, branch cuts, frac¬ 
tional powers, discontinuities in the function or in any of 
its derivatives etc. in the complex plane. Important in 
this context is the convergence domain of a Chebyshev 
series. It is given by the interior of an ellipse whose foci 
lie at X = ±1. This is in contrast to a Taylor series (or 
more generally a Laurent series), whose domain of con¬ 
vergence is a disc around the expansion point. 

We shall close this mini-review by recalling the Cheby¬ 
shev truncation theorem, which gives an upper bound for 
the error made in truncating a Chebyshev series. The 
error is given by the sum of the absolute value of the 
neglected coefficients. A useful rule of thumb is that this 
error is of the order of the last retained coefficient for 
exponential convergence, and of the order of the num¬ 
ber of coefficients retained times the last coefficient for 
algebraic convergence. 

III. IMPLEMENTATION 

In this section, we want to give more details on how to 
apply pseudo-spectral methods by considering a generic 


f{x) = 00 - 1 - xbi - 62 . 
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ODE of one function in one variable. Without loss of 
generality we restrict to the domain x G R+. If the 
domain is R, then we apply the subsequent ideas to both 
R+ and R_. 

Let £ be a (not necessarily linear) integro-differential 
operator and consider the problem 

£[f(x)]= 0 . ( 8 ) 

We want to decompose the function f(x) into a series 
of Chebyshev polynomials. To gain maximum efficiency, 
first the domain of f(x) is decomposed into M parts. We 
will restrict the discussion here to Af = 2, that is we 
decompose the domain into [0,a;o] and [a;o,oo]. In the 
first domain, the function is interpolated via a standard 
Chebyshev series, whereas in the second part, a rational 
Chebyshev series is used. Thus 

Z; -1), x<Xo, 

N. (9) 

foo(x) Y. riRi{x - xo), X > Xo, 

i=0 

where fao{x) is the leading term of the asymptotic be¬ 
haviour of f{x) for X -G- oo, which can be easily deter¬ 
mined analytically in many cases. For FRG equations, 
the asymptotics is typically determined by dimensional 
scaling properties. This ansatz can be inserted into ( 8 ). 
To solve for the coefficients, it is useful to apply the collo¬ 
cation method. It consists of evaluating the equation on 
a certain set of collocation points and solving the result¬ 
ing algebraic set of equations, e.g. by a Newton-Raphson 
method. The key to high accuracy is the choice of these 
collocation points. It turns out that the best choice is 
to employ either the nodes or the extrema of Chebyshev 
polynomials [ 86 ] . Additionally, one must match the func¬ 
tion value as well as derivatives of both expansions at the 
intermediate point Xq to achieve smoothness. If the dif¬ 
ferential equation is of order p, then p— 1 derivatives have 
to be matched. 

It should be noted that in the above decomposition, 
one has two free parameters: the matching point xq as 
well as L, which encodes the compactification of the semi¬ 
infinite domain. There is no intrinsic rule how to choose 
them. We found it reasonable to choose Xq large enough 
to include the essential physics, e.g. the vacuum expec¬ 
tation value in a scalar field theory. The influence of L 
is usually small in this case, as then the rational Cheby¬ 
shev expansion essentially only interpolates the asymp¬ 
totic behaviour. Either way, if the expansion converges, 
it will converge for any choice of the parameters. 

The above ideas can be generalised in two ways. 
Firstly, one can trivially apply these methods to a system 
of functions, i.e. a coupled system of ODEs. An exam¬ 
ple of this will be given later, and can also be found in 
[57]. Secondly, via a tensor product, the generalisation to 
PDFs is possible. This will be addressed in future work. 


IV. THE FUNCTIONAL RENORMALISATION 
GROUP 


A very efficient means to store the full quantum in¬ 
formation of a quantum field theory is the so-called ef¬ 
fective action T, which is defined as the Legendre trans¬ 
form of the Schwinger functional. There are numerous 
possibilities to compute T, one is given by the FRG. In¬ 
stead of r, the effective average action T^ is considered, 
which smoothly interpolates between a microscopic the¬ 
ory Tfe^A = 5'cb where A is an ultraviolet cut-off, and the 
full quantum theory Tfe^g = T- Following Wilson’s idea, 
quantum fluctuations at momentum scale p^ ~ k'^ are 
successively integrated out during this evolution. This 
process is described by an exact FRG equation, the Wet- 
terich equation [ 1 ], 


kdkTk = -STr 


\kdkRk) 


( 10 ) 


( 2 ) 

where T), denotes the second functional derivative of 
Tfc with respect to the fields and the super-trace STr 
stands for a summation over discrete indices, integration 
over continuous indices and an additional minus sign for 
Grassmann-valued fields, i.e. fermions. The functional 
Rk is a regulator, which ensures both infrared as well as 
ultraviolet Hniteness. Detailed information on the FRG 
can be found in e.g. [3, 27, 28, 62]. 

In most cases, the exact functional integro-differential 
equation ( 10 ) can only be solved by choosing a certain 
truncation for the effective average action. A class of 
common truncations is the derivative expansion, which 
takes derivative interactions up to a given order into 
account. For many cases, such a systematic expansion 
yields results comparable to those obtained by other 
methods, e.g. lattice studies. 


V. 0(1) MODEL NEAR CRITICALITY 

This section is devoted to a detailed study of various 
properties of the 0(1) model. Our ansatz for the effective 
average action reads 

LfcM = y d^'x |izfe(cr^)(i9^cr)^-t t/fe(cr^)| , (11) 

which contains the effective potential Uk{<J^) and a wave 
function renormalisation Zk{cr^). The operators are cho¬ 
sen such that the symmetry of the scalar field cr is pre¬ 
served. In first order derivative expansion, also called lo¬ 
cal potential approximation (LPA), one neglects the run¬ 
ning and the field dependence of the wave function renor¬ 
malisation, Zkicr'^) = 1. By contrast, within next-to- 
leading order in the derivative expansion (NLO), the full 
flow of the wave function renormalisation is taken into 
account. As a compromise between LPA and NLO, one 
often considers a field-independent but scale-dependent 
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FIG. 1: Derivative of the effective potential of the Wilson- 
Fisher fixed point in LPA and LPA’. 

wave function renormalisation Z}. = ^^((Tg), usually 
called LPA’. Here, cto is typically chosen to be the vac¬ 
uum expectation value of the scalar field. 

Here, we first study the Wilson-Fisher fixed point in 
three dimensions in LPA and LPA’. Although there ex¬ 
ists a full analytic solution for the large-case, we will 
confine ourselves to the case of = 1. A large-Ai study 
is given below for the Gross-Neveu model in section VI. 
Subsequently, we look for multi-critical fixed points in 
dimensions 2 < d < 3 [11]. 

A. Wilson-Fisher fixed point in LPA and LPA’ 

First, let us study the well-known Wilson-Fisher fixed 
point of the 0(1) model in three dimensions. In arbitrary 
dimension d, the fixed point equation is given by [2] 

(-2 -h T])u'{p) + {d-2 + ■n)pu''{p) 

4^’d A d ^ 3 m"(p)- h2pu'"(p) _ 

d 'v d + 2) {l + u'{p) + 2pu"[p)Y ’ ^ ’ 

where u{p) = is the dimensionless effective po¬ 

tential as a function of the dimensionless invariant p = 
= 2'^+i7r‘^/2r(d/2) and p = -kdklnZk 
is the anomalous dimension. The latter is given by [87] 

^ l&Vd pQu''{pof 

^ d (1 + 2pou"(po))2 ’ 

po being the vacuum expectation value (vev). In these 
equations the optimised regulator is employed [63]. For 
aspects of optimisation, see also [64]. As (12) does not 
depend on the potential itself, all our calculations will 
involve its derivative instead, i.e. f{x) u'(p). 

In the following, we will compare the solution to (12) 
in LPA, i.e. with p = 0, and in LPA’, where we include 
the anomalous dimension, but no field-dependent wave- 
function renormalisation. As numerical parameters, xq = 


3/10 and L = 1 were chosen, and we used floatl28 
which gives twice as many figures as the conventional 
double data type. From (12), one infers the asymptotic 
behaviour u'^{p) = 

Figure 1 displays the derivative of the effective poten¬ 
tial. One can see that the inclusion of the anomalous 
dimension has a quantitative influence for intermediate 
values of p. For the vev, we get 

pLPA ^ 0.030647942408697774953, 

pLPA’ ^ 0.030592776234779436405. (14) 

In LPA’, we found the anomalous dimension to be 

p = 0.044272337370315035214. (15) 

It may seem ridiculous to present that many figures, 
but they illustrate the power of pseudo-spectral methods. 
Our viewpoint here is to solve a truncated problem (nu¬ 
merically) exactly, and any numbers given here are to be 
understood as the solution to the truncated problem. 

We find that our values match very well with earlier 
results, e.g. given in [3, 5, 9], and [65] where fixed point 
quantities were calculated via full potential flows. In par¬ 
ticular, we can reproduce all digits of the high precision 
results given in [6], where an LPA truncation was used 
(notice the rescaling of p by 87r^ compared to our con¬ 
ventions). 

In comparison to results obtained from High- 
Temperature expansions and Monte-Carlo simulations 
(p = 0.036) [66], our value has a systematic error which 
is to be expected in the LPA’ truncation. Within the 
FRG, an accuracy competitive with other methods has 
been reached using the BMW approximation technique 
[8, 24]. 

Figure 2 shows the coefficients of the expansion of u'{p) 

2-r) 

in a Chebyshev series and of u'{p)/p^+-^ in a rational 
Chebyshev series. In LPA, one perfectly sees exponen¬ 
tial convergence in both the Chebyshev as well as the 
rational Chebyshev series. On the other hand, as soon as 
we include the anomalous dimension, we find asymptot¬ 
ically only algebraic convergence in the rational Cheby¬ 
shev case. This behaviour is indeed expected by the 
asymptotic behaviour of the potential, as it rises with 
a fractional power. Furthermore, one can also see that 
this problem is irrelevant for all practical purposes, as 
the algebraic convergence only sets in at about 10“^®, up 
to that point one still observes exponential convergence. 
This emphasises the fact that any statement about con¬ 
vergence is really an asymptotic one, and one cannot pre¬ 
dict where this behaviour sets in. As a final comment on 
this, note also the number of coefficients needed to gain 
a certain accuracy: in case of exponential convergence, 
one needs very few coefficients to get an adequate result, 
but as soon as there are singularities of any kind, one 
needs a large number of coefficients to further increase 
the accuracy. 
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LPA LPA' 



FIG. 2: Decay of coefficients of Chebyshev and rational Chebyshev expansion of the derivative of the Wilson-Fisher fixed point 
potential, LPA on the left and LPA’ on the right. Notice the algebraic decay in the latter case in the rational Chebyshev region, 
which needs a factor of 10 as many coefficients as in the case of exponential convergence in LPA to achieve the same order of 
accuracy. 


Another point of interest is that no additional condi¬ 
tion has to be imposed in solving (12), in particular no 
boundary condition or the like. Only the asymptotic be- 

2 — rj 

haviour u'{p) oc and a sufficiently good initial guess 
is needed. This may seem unexpected for a differential 
equation of second order, but can be understood along 
the lines of [67]. Indeed, analysing the situation, one can 
see that at p = 0, the order of the differential equation 
decreases by one, which fixes one condition, and the same 
is true at p = oo. 

As a further case for the method, let us expand the 
solution in LPA into a Taylor series around vanishing 
field as well as a Laurent series around p = oo and 
compare whether the relations between coefficients ob¬ 
tained by plugging in such an ansatz into the hxed 
point equation are satisfied. Around vanishing field, with 
w*(p) = one obtains the well-known relations (see 

e.g. [23]) 

Oi = —47r^ao(l -|- Oq)^ , 

12 

a 2 = —'7r^ao(l + ao)^(l + 13ao), 

5 

OOQ 

03 = 7r®Oo(l -b ao)‘‘(l -b 7ao), 

32 

04 = —7r®OQ(l -b ao)®(2 -b oo(121 -b 623oo)), (16) 


etc. Inserting our solution, one finds that the absolute 
error in these coefficients are (< 10“^°, 2 x 10“^^, 2 x 
10“^®, 7 X 10“^®). For the expansion around infinity, one 
finds that 


<(p) = V 


1 

75bl7r2p3 


0{p-^). 


(17) 


Expanding our solution, the coefficients of 


p^, p^, p~^, p~^ (which should vanish in the exact 
solution) are (—4 x 10“^^, 3 x 10“^'^, —8 x 10“^^, 10“^®), 
and the relation between the leading and the first 
sub-leading coefficient is fulfilled to an absolute accuracy 
of 10“^^. For completeness, let us give the values of oq 
and A both in LPA and LPA’: 

= -0.18606424947031443565, 

= -0.16574071049155738982, 

^LPA ^ 84.182303273336100651, 

^ 50.323366981670544177. (18) 

These results match with [23] and [68] where local ex¬ 
pansions and the shooting method were employed. This 
underlines that we can trust the global solution and that 
we can relate to earlier results. 

Let us now turn our attention to the critical expo¬ 
nents of the Wilson-Fisher fixed point. They are defined 
as minus the eigenvalues of the linearisation of the per¬ 
turbed fixed point equation. Again, a global approach to 
the solution of the perturbed equation is used. Figure 3 
shows the eigenfunctions corresponding to the five high¬ 
est eigenvalues, where the anomalous dimension has been 
taken into account. As for the potential itself, any pre¬ 
cision can be achieved in the eigenfunctions and critical 
exponents. The critical exponents match with earlier re¬ 
sults, e.g. given in [.3]. Having said that, let us emphasise 
again that the largest error arises from the systematic er¬ 
rors of the derivative expansion to order LPA/LPA’. If 
we compare with Monte-Carlo results [66], we find a de¬ 
viation of about 2.5% for the first and 27% for the second 
critical exponent. Especially the error of the second crit¬ 
ical exponents is to be expected from the low order of the 
derivative expansion used, see [5]. 
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e= 1.548 

4 - e=-0.6127 

. e=-3.073 

- e=-5.738 



FIG. 3: Eigenperturbations of the Wilson-Fisher fixed point, 
normalised to one at p = 0. 

As shown in this example, the error is dominated by 
truncating the effective average action and not by numer¬ 
ical errors. For this reason, from now on we will only give 
a few relevant digits, bearing in mind that in principle 
we could calculate as many digits as needed. 


B. Multi-critical fixed points for 2 < d < 3 

It is worthwhile to have a closer look at fractional di¬ 
mensions 2 < d < 3. The fixed point structure is get¬ 
ting richer for decreasing dimension and, therefore, it is 
interesting to investigate the interpolation between the 
two fixed points in d = 3, the Gaussian and the Wilson- 
Fisher fixed point, and the infinite number of fixed points 
in d = 2. In [11, 69] the existence and properties of 
multi-critical fixed points in dependence on d and N are 
investigated. These results are used to give an RG proof 
of the Mermin-Wagner-Hohenberg theorem [70-73]. 

As a test case and as has been done in [9] we restrict 
ourselves to the Ising universality class = 1 here. We 
emphasise that the following investigations can straight¬ 
forwardly be applied to arbitrary flavor numbers N if 
the set of fixed point solutions is still discrete. In [9], a 
sequence of critical dimensions dc,i where the next multi- 
critical fixed point potential Ui (p) emerges was suggested. 
These are those dimensions where new operators p’’ be¬ 
come relevant for d < d^^i- Concentrating in the fol¬ 
lowing on d = 2.4 as an example, we find three more 
multi-critical fixed points FPjg{ 3 ^ 4 ^ 5 } besides the Wilson- 
Fisher (FP 2 ) and the Gaussian (FPi) fixed point. The 
index i > 2 counts the minima of the corresponding fixed 
point potential counted in the dimensionless scalar field. 
Our results around d = 2.4 confirm the predicted value 

— 5 ■ 

For our calculations we have employed (12) and (13) 
within the LPA’ truncation. The anomalous dimension 
is again evaluated at the global minimum of the potential 
which is in the following cases the outermost minimum. 
In Figure 4 the first derivative of the multi-critical fixed 



WF-FP 


V 

relev. exp. 

irrelev. exp. 

0.1390 

1.1441 

-0.7919 

-3.1129 

-5.6370 

multi-critical FP 

i—3 

n 

relev. exp. 

irrelev. exp. 

0.01598 

1.9629 

-0.5108 


0.8416 

-2.0698 

-3.8140 

multi-critical FP 

i=4 

V 

relev. exp. 

irrelev. exp. 

0.001753 

1.9969 

-0.3138 


1.4615 

0.6726 

-1.3968 

multi-critical FP 

i—5 

V 

relev. exp. 

irrelev. exp. 

8.2715 X lO"'" 

1.9999 

-0.1655 


1.5973 

1.1243 

0.5414 

-0.9297 


TABLE I: Anomalous dimensions and highest critical expo¬ 
nents of all scaling solutions in d = 2.4. 

point potentials are shown. As the values of the anoma¬ 
lous dimension of the multi-critical fixed points FPi >3 
are small compared to the one of the Wilson-Fisher fixed 
point, the convergence of their coefficients is exponential 
within the used precision of float 128. Therefore, the 
deviation from the exact solution can be estimated to be 
below I0“^°. 

In Table I the anomalous dimensions and the largest 
critical exponents calculated by pseudo-spectral meth¬ 
ods are given. Our results are in good agreement with 
[9, 11, 69]. Additionally, the results for the Wilson-Fisher 
fixed point in d = 2.4 can be related to earlier works [74- 
76], where the £-expansion and lattice simulations were 
applied. As already noted in [11], the highest relevant 
critical exponent for the fixed points i >3 is close to the 
mean field value 2 at the corresponding critical dimen¬ 
sion. The other relevant exponents are smaller. 

The sequence of critical dimensions predicts that a new 
fixed point potential with six minima (regarded as func¬ 
tion of the dimensionless scalar field) emerges exactly at 
d = 2.4. As the Wilson-Fisher fixed point probably does 
not exist in d = 4 but exists in all dimensions 2 < d < 4 
we find this fixed point for all dimensions d < 2.4. In fact, 
we are able to determine a global solution for d = 2.399 
where the non-asymptotic behaviour is realised on very 
smalls scales \u'{p < po)\ ^ 10“® and rj = 2.3446 x 10“^°. 


VI. GROSS-NEVEU-MODEL IN d = 3 

In this section we extend our studies to the par¬ 
tially bosonised Gross-Neveu model in d = 3 dimen¬ 
sions. Loosely speaking, this is a generalisation of the 
0(1) model including N fermionic degrees of freedom. It 
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FIG. 4: First derivative of multi-critical fixed point potentials exhibiting two minima regarded as function of the dimensionless 
scalar field (blue, Wilson-Fisher potential), three minima (yellow), four minima (red), five minima (green). The small insets 
depict the global behaviour of the solutions. 


has applications in condensed matter physics and serves 
as a toy model for asymptotic safety scenarios. A detailed 
analysis can be found in [15]. 

The ansatz for the effective action in an LPA’-type 
truncation reads 

Tk['tp,'tp,cr] = / + ihka)ip 

^ (19) 

+ + [/(cr^)| . 

The bosonic (Za-^k) and fermionic (Z^^k) wave function 
renormalisations, and the Yukawa coupling hk, which 
transfers the interaction between bosons and fermions, 
are assumed to be scale-dependent but field-independent. 
In the following studies, we again employ a Litim-type 
cutoff [63]. 


A. Large-A analysis 

The large-A approximation is a good test case because 
the fixed point equations can be solved analytically. In¬ 
terestingly the scalar anomalous dimension does not van¬ 
ish in contrast to the one of the Wilson-Fisher fixed point. 


Even the fixed point potential looks very different. The 
fixed point equations in the large-A limit are given by 
[15] 


0 =(-2 -b ri^)u'{p) + id-2 + ri^)u"{p)p 


M^Vd A _ pp 


d+lj (l-b2/i2p)2 


Q ={^d - A + 2pp + p^)h? 


d'y'^d I ^ Vp 


Pa 

d 


4^ 2d-A 


Pp = 0 , 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


again denoted in dimensionless quantities u{p) = 
k-‘^U{a‘^)/N,_ p = Z<,,fcfc2-^aV(2A) and = 

Z~l.Zz\k'^~'^h'^N^ and where d^ stands for the dimen¬ 
sion of the Dirac algebra. Note that for the large-A limit, 
an appropriate rescaling has been taken into account. 

In this approximation, we encounter a first order sys¬ 
tem. The bosonic anomalous dimension can be read off 
from (21) to be ?7o. = 1 exactly (as long as h ^ 0). We 
can reproduce this result to all digits of float 128 which 
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FIG. 5: Decay of coefficients of Chebyshev and rational 
Chebyshev expansion in the large-limit of the Gross-Neveu 
model. 



gives an accuracy of about 10 The exact fixed point 
value of the Yukawa coupling reads 


h 


2 

* 


d \ (d-4)(d-2) 

d-yVd ) (8 — 6d) 


(24) 


and can be confirmed up to 10“^^ as well. The fixed 
point potential is given by the Gaussian hypergeometric 
function [15], 


u^{p) = 


4(8-6d- h d^) 
3d-^4 


px 


/ 1 2 — d d 8 — 6d-|-d^\ 

\1 —d’ ’ 1 — d’d-yUd 3d — 4 / 


(25) 


The absolute difference between the analytic solution 
and our numerical one can be estimated to be smaller 
than 3 x 10“^^ for large p. For finite p it is even smaller. 
This is due to the Gaussian grid which only has points 
at finite p. Thus the asymptotic prefactor is only tuned 
regarding finite field values and, therefore, has a larger 
error of about 3 x 10“^^. For this calculation we have 
used xq = 3/10 and L = 2. The decay of the coeffi¬ 
cients can be seen in Figure 5. The Chebyshev expansion 
shows exponential convergence. By contrast, the ratio¬ 
nal Chebyshev coefficients decrease exponentially at first, 
but only up to a certain number of coefficients, the ac¬ 
tual convergence rate is algebraic. This is to be expected 
due to the asymptotic behaviour oc ^/x. The behaviour 
of the last coefficients shows a truncation effect which is 
not a numerical effect. If we calculate the spectral co¬ 
efficients from the analytic solution, we actually obtain 
a good agreement with the numerically calculated ones. 
Ignoring the last coefficients affected by the truncation 
we read off ^ 10“^® for the lowest coefficient. The rule 
of thumb that the error can be estimated by iV • cat is 
in very good agreement with the maximal deviation of 
about 3 X 10“^^ from the exact solution. 



20 40 60 80 100 


n 

FIG. 6: Decay of rational Ghebyshev coefficients in depen¬ 
dence on L in the large-Y limit of the Gross-Neveu model. 


Let us have a closer look at the choice of the two pa¬ 
rameters Xq and L. They have quite some influence on 
the convergence behaviour of our coefficients. We observe 
that if the matching point Xg is chosen to be smaller, then 
the decrease of the Chebyshev coefficients is exponential 
as well, but faster. This also holds for the rational coeffi¬ 
cients if the matching point is increased. However, lower¬ 
ing Xq the algebraic convergence sets in earlier, enlarging 
Xq the algebraic convergence sets in later. It is remark¬ 
able that the gain of accuracy is inappreciable when tak¬ 
ing enough coefficients into account because if algebraic 
convergence has set in, the coefficients do not change sig¬ 
nificantly. This is different for the parameter L. Figure 6 
shows that one can gain orders of magnitudes of accuracy 
if L is increased. The decrease of the first coefficients is 
more slowly, but the algebraic convergence sets in later. 
This short analysis already makes clear that the choice of 
optimised parameters can strongly depend on the max¬ 
imal number of coefficients that one takes into account. 
These observations from one specific example may give 
an indication for other calculations as well. 


B. Finite N analysis 

For finite IV, the fixed point potential shows some in¬ 
teresting behaviour. In [15] it is shown that the fixed 
point potential lies in the symmetric regime for all iV > 2. 
For that reason, the potential was expanded polynomi- 
ally. A study of the convergence radius indicates the 
reliability of these results. Unfortunately, the conver¬ 
gence is less clear for smaller N such that a global solu¬ 
tion is required. For small JV the fixed point potential 
moves from the symmetric regime to the spontaneously 
symmetry-broken regime. In [77] a fixed point potential 
for N = ^ (corresponding to one Dirac fermion in the 
irreducible representation) was found in the symmetry- 
broken regime. In [15] it was assumed that the non- 
Gaussian Gross-Neveu fixed point interpolates between 
the large-IV fixed point and the Wilson-Fisher fixed point 
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FIG. 7: Fixed point potential (npper left panel), relevant exponent (npper right panel), Ynkawa-coupling (lower left panel) and 
anomalous dimensions (lower right panel) for the Gross-Neveu model with 0.3 < Y < 12. For increasing flavor number u(,(0) 
increases as well. The dashed line depicts the flavor number Nt of the transition. 


in the iV —>■ 0 limit. 

The hxed point equations for the Gross-Neveu model 
read, for general values of N [78], 


0 =(-2 -I- r]a)u'{p) + {d-2 + ■q„)u"{p)p 


(l- ^ 3u"(p) -I- 2pu"'{p) 

d V d + 2) {I + u'{p) + 2pu"{p)Y 

M^Vd ( _ Tip \ 

d \ d+iy (l + 2h2p)2 ’ 


(26) 


0 - A + 2r]^ + ri^)h'^ 




1 _ 

^ d+1 


+ 


1 _ 

^ d+2 


l + 2poh'^ l + u'{po)+ 2pou''{po) 


(1 + 2po^^) + u'{po) + 2pou"{po)) ^ 


{l + 2poh'^) ^{1 + u'{po) +‘ipou'^po)) 


, 6 „ [r, ^ d+1 


-h^Po 2 


1 _ V,T 

^ d+2 


l + 2poh'^ l + u'{po)+ 2pou''{po) 


(1 + 2poh^)-^{l + u'ipo) + 2pow"(po))-', (27) 

=S^h^Nx 


1 I 

4 2d-4 


-h 


d-2 


IP =- 


1 — 2poh^ 

(1 + 2poh2)4 ' (l + 2poh2)2 ' (l + 2poh2)3^ 

(28) 
, (29) 


8vd (3m"(po) -h 2pou'"(po))2 
d (1 -h u'{po) + 2pou"{po))'^ ’ 


d (1 -I- 2poh2)(l -I- u'(po) + 2poM"(po))^ ’ 


- (48 pom"(po) + 32 pom'"(po)) x 

/ 1_1_!Z£- \ 

I _ d-\-\ I 2_ d-\-‘2 _ \ 

yl + 2poh2 1 -I- m'(po) + 2pou"{po) j 


with u{p) = k ‘^U{a‘^), p = (2 and h? = 

Z~\Z^\k'^~'^h'^■ The asymptotic behaviour of the po- 

2 — ri(j 

tential is given by u'{p) oc pd~ 2 +'n^ , 

We have calculated the global solution to these hxed 
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FIG. 8: Global fixed point solution for the scalar-tensor model in three dimensions. The scalar potential closely resembles the 
Wilson-Fisher fixed point potential on flat background. The non-minimal coupling is strictly positive, admitting a positive 
Newton constant. 


point equations and the relevant critical exponent in 
three dimensions for N lying between 0.3 and 12, see 
Figure 7. We obtain a very good agreement with [15] for 
N > 2 where a polynomial approximation is employed. 
Even the relevant exponent matches in the first four rele¬ 
vant digits. It is worth mentioning that this good agree¬ 
ment is only obtained by taking high orders in the poly¬ 
nomial truncation into account [15], especially for small 
N. Our results for N >2 are also compatible with other 
methods such as 1/iV-expansions [79, 80], and Monte- 
Carlo simulations [81, 82]. In fact, systematic truncation 
errors appear to be smaller for the Gross-Neveu model in 
comparison with the 0(1) model. The overall consensus 
among the non-perturbative methods is very satisfactory. 

Let us now concentrate on the small-A'^ regime. The 
transition from the symmetric to the symmetry-broken 
regime can be determined to be at Nt ~ 0.5766. As a 
new result, we observe that the Gross-Neveu fixed point 
does not approach the Wilson-Fisher fixed point for small 
N. This can be seen from the behaviour of the Yukawa 
coupling and the anomalous dimensions on the one hand 
and the relevant exponent on the other hand. In partic¬ 
ular, the behaviour of h* suggests that the Gross-Neveu 
fixed point moves to infinity in theory space for N ^ 0. 

It is instructive to compare our results for = 1/2 
with those of [77] where also a full potential flow has 
been studied (note that our convention of Af = 1/2 cor¬ 
responds to Af = 1 in [77]; for aspects of criticality 
see [20]). In [77] the fluctuation terms oc pou''{po)h'^, 
oc pQu"'{po)h'^ and cx pgh^ have been missed in the 
derivation of the flow equation, see the discussion in [56] . 
If we artificially switch off these terms, the vacuum expec¬ 
tation value and the critical exponent of our calculation 
are in good agreement with those of [77]. On the con¬ 
trary, including these terms, even the first relevant digit 
changes. For N = 1/2, we obtain v = l/0i = 0.4836, 
77 cr = 0.3227, 77 ^ = 0.1204. In conclusion it is remarkable 
that our approach is able to find a global solution in a 
regime where a polynomial truncation is not reliable. 


VII. SCALAR-TENSOR GRAVITY 

As a final example, let us consider a model which cou¬ 
ples a scalar field non-minimally to gravity [50, 83, 84]. 
The ansatz for the effective average action is given by 

^k[(J,g]= J (vk{cr) - Fk{a)R+^g^"''df,ad^a^ , 

(30) 

where g denotes the determinant of the metric < 7 ^^ and R 
is the Ricci scalar. It serves as an effective model for the 
cosmological evolution of the Universe. Here, we want to 
use the flow equations of [50], in which an exponential 
split was used to quantise the gravitational fluctuations. 
The explicit flows of the dimensionless variants v{p) and 
f{p) of the functions V(a) and F{a) can be looked up in 
[50]. 

We have found a non-trivial solution to the fixed point 
equations, which is shown in Figure 8 . Some remarks are 
in order. First, the effective potential of the scalar field 
closely resembles the Wilson-Fisher fixed point potential 
on flat background. Second, note that /»(p) > 0 for all 
p > 0, which implies that the (analogue of the) Newton 
constant is positive, and we are indeed in the physical 
regime. Third, the minimum of the potential v{p) lies at 
p = 0.05004. Lastly, we observe exponential convergence 
of the expansion. Thus, the solution can in principle be 
computed to arbitrary precision. 

Let us now turn our attention to the critical exponents 
of the fixed point. The six leading critical exponents are 

01 = 3, 02 = 1.9134, 

03 = 1.1798 , 04 = 0.6679, (31) 

05 = -0.2812, 06 = -1.217. 

Notice that the first exponent is present in any fixed point 
solution. It is related to what one would call cosmologi¬ 
cal constant, being the most relevant perturbation. The 
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FIG. 9: Eigenfunctions of the fixed point solution to the three-dimensional scalar-tensor model. They are normalised such that 
5u(0) = 1. From relevant to irrelevant critical exponents: thin blue, thick orange, dotted green, dashed red, dot-dashed violet. 


eigenfunctions for all of them, except for the first one, 
are given in Figure 9. For the first critical exponent 0i, 
the eigenfunction is {Sv = l,Sf = 0). 

Let us compare these findings with those of [50] . There, 
the fixed point equations were first analysed in a one- 
loop approximation. Using a local field expansion, two 
non-trivial solutions have been identified: the conformal 
solution V = l/(187r^),/ = 37/(727r^) -|-(t^/ 4, possessing 
four relevant directions, and a Wilson-Fisher-like solution 
with three relevant directions. There are indications that 
the latter fixed point runs into a singularity due to / 
having a zero. On the basis of the full equations, no 
conformal solution with the same simple structure as in 
one-loop approximation was found. No statement was 
given on the fate of the Wilson-Fisher-like solution. 

Our results are the following: the critical exponents of 
our global solution are very close to the ones of the con¬ 
formal solution in the one-loop approximation. Also, the 
general form of the non-minimal coupling f{p) is quali¬ 
tatively the same. On the other hand, the form of the 
scalar potential v{p) is Wilson-Fisher-like. 

We further checked deviations from three dimensions. 
For any dimension between three and four, we find a sim¬ 
ilar fixed point solution. Though, at d = 4 exactly, this 
solution seems not to be present anymore, reminiscent of 
the situation in a purely scalar theory. We take this as 
an indication that the solution found here might be the 
generalisation of the Wilson-Fisher fixed point. 

Finally, let us comment on why solutions, where 
the non-minimal coupling becomes negative, f{p) < 0, 
should be taken with care. In the derivation of the flow 
equation, the regulator was spectrally adjusted. In par¬ 
ticular, for the tensor fluctuations, the regulator was cho¬ 
sen to be proportional to f{p). As soon as this function 
crosses zero, the tensor modes are not properly regu¬ 
larised anymore. In the flow equations, this is reflected 
by terms proportional to l//(p), which can only be com¬ 
pensated by divergences in derivatives of the potential. 
This might explain why we could not find further global 


solutions. 

Summarising, we find a globally well-defined solution 
in scalar-tensor gravity. However, the inclusion of matter, 
in particular scalar fields, induces higher order curvature 
terms. Thus, one should check that the global solution 
survives. On the other hand, for cosmological applica¬ 
tions, an Einstein-Hilbert type of truncation as discussed 
here is expected to be well-suited as an effective model. 


VIII. CONCLUSIONS 

In this work, we presented a method to solve ODEs 
globally. Pseudo-spectral techniques are not new and al¬ 
ready have been applied to various problems in physics. 
A lot of rigorous results on spectral or pseudo-spectral 
methods are available [49]. However, to our knowledge 
the expansion in rational Chebyshev polynomials is not 
very well established in quantum field theory calcula¬ 
tions, although it allows to investigate the question of 
global existence of solutions. This is a very important 
question since the non-linear ODEs encountered in ERG 
studies can have many more or less stable local solutions 
which are not easy to distinguish from global ones if only 
local information is accessible. For instance, the physical 
criteria of polynomial boundedness and self-similarity are 
difficult to impose locally [85]. The method presented in 
this paper offers a comparatively easy way to find global 
solutions of ODEs that circumvents such pitfalls. 

We applied this method to various models. The first 
test case was the very well known 0(1) model in three 
dimensions which we considered in both the LPA and 
LPA’ truncation. There are numerous works on expan¬ 
sions for small and large fields and results gained via 
the shooting method which give a good impression of the 
global behaviour of the Wilson-Fisher fixed point. We re¬ 
ported on the difference between LPA and LPA’ trunca¬ 
tions taking the global behaviour of the potential into ac¬ 
count. Although the anomalous dimension is very small, 






12 


the asymptotic behaviour, especially with regard to the 
prefactor, changes significantly. Besides the fixed point 
potential itself, we calculated the eigenfunctions globally 
and determined the critical exponents. For all quantities 
we obtained good agreement with already known results 
calculated with other methods. As far as numerical ac¬ 
curacy is concerned, our method outperforms previous 
results by many orders of magnitude, while being very 
stable, fast and lightweight. 

Subsequently, we extended our study to fractional di¬ 
mensions, taking d = 2.4 as a representative. We found 
all multi-critical fixed point potentials predicted in [9] 
and could, moreover, determine their global behaviour. 
We were able to see the next higher critical fixed point 
emerging at d < 2.4 which demonstrates that our numer¬ 
ical method is highly accurate and stable. All physical 
quantities, the anomalous dimension and critical expo¬ 
nents, again match with earlier results. 

As a next system we considered the Gross-Neveu 
model. On the one hand, the large-limit provides an 
easily accessible analytic solution. On the other hand, 
the small-iV limit is not easily accessible by use of com¬ 
mon local expansions and offers, therefore, the possibil¬ 
ity to demonstrate the advantages of our global method. 
In the large-A^ case we obtained a conclusive agreement 
with analytic results. For finite N, our results agree very 
well with other data, including 1/iV-expansions and lat¬ 
tice methods. We were able to fix the transition fla¬ 
vor number to be Nt ~ 0.5766 and observe how the 
fixed point potential goes over from the symmetric to the 
symmetry-broken regime. We found that the fixed point 
Yukawa coupling grows large for N —>■ 0. The anoma¬ 
lous dimensions, in particular the fermionic one, take on 
finite values. This suggests that the Gross-Neveu fixed 
point does not merge with the Wilson-Fisher fixed point 
in the limit Y —>■ 0 contrary to what has been anticipated 
in [15]. Additionally, we saw that all fluctuation terms 
in the Yukawa fixed point equation which occur in the 
symmetry-broken regime have a significant influence on 


physical quantities, such as critical exponents. Compar¬ 
ing to [77] where some fluctuation terms were missed we 
determine the deviation to be up to 30%. 

We finally discussed a scalar-tensor gravity model in 
d = 3. This model is supposed to have a gravitationally 
dressed Wilson-Fisher fixed point. We found a global 
solution implying a positive Newton constant with four 
relevant critical exponents. No other non-trivial global 
solution was found, though it may be difficult to prove 
that our solution is unique. We could recover the solution 
in all dimensions smaller than four, but not in d = 4. 
This is taken as an indication that the solution found 
might be the generalisation of the Wilson-Fisher fixed 
point to curved space. 

We emphasise that the models treated in this pa¬ 
per only represent a small set of examples for a large 
amount of possible applications of pseudo-spectral meth¬ 
ods. These methods are straightforwardly extendable to 
include more than one variable which is needed if addi¬ 
tional invariants are taken into account. Full potential 
flows requiring the technique to be extended to PDFs 
are underway. 
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